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Abstract 

We revisit the observational constraints on compact (closed) hyperbolic(CH) models 
from cosmic microwave background(CMB). We carry out Bayesian analyses for CH 
models with volume comparable to the cube of the present curvature radius using 
the COBE-DMR data and show that a slight suppression in the large-angle temper- 
ature correlations owing to the non-trivial topology explains rather naturally the 
observed anomalously low quadrupole which is incompatible with the prediction of 
the standard infinite Friedmann- Robertson- Walker models. While most of positions 
and orientations are ruled out, the likelihoods of CH models are found to be much 
better than those of infinite counterparts for some specific positions and orientations 
of the observer, leading to less stringent constraints on the volume of the manifolds. 
Even if the spatial geometry is nearly flat as Qtot = 0.9 — 0.95, suppression of the 
angular power on large angular scales is still prominent for CH models with volume 
much less than the cube of the present curvature radius if the cosmological constant 
is dominant at present. 



1 Introduction 

In the framework of modern cosmology, one often takes it for granted that the spatial 
geometry of the universe with finite volume is limited to that of a 3-sphere. However, 
if we drop the assumption of simply connectivity of the spatial geometry, then flat and 
hyperbolic (constantly negatively curved) spaces can be closed as well. Therefore non- 
positive curvature of the spatial geometry does not necessarily implies the infinite space. 
For instance, the simplest example of a constantly curved finite-volume space with non- 
trivial topology is a fiat 3-torus which can be obtained by identifying the opposite faces of 
a cube by translations in the Euclidean 3-space. Construction of compact (closed) hyper- 
bolic(CH) spaces is somewhat complicated but one can systematically create a countably 



1 



infinite number of topologically distinct classes of CH spaces by performing topological 
surgeries on a certain space. In order to avoid confusion, it might be better to call a 
cosmological model whose spatial geometry is described by a connected and simply con- 
nected constantly negatively curved space as an "infinite hyperbolic" model rather than 
an "open" model which has been commonly used in the literature of astrophysics. 

Since 1993, a number of articles concerned with constraints on the topology of flat 
models with no cosmological constant using the COBE-DMR data have been published. 
[05 il; H) @) Hi The large-angle temperature fluctuations in the cosmic microwave back- 
ground(CMB) discovered by the COBE satellite constrain the topological identification 
scale L Q larger than 0.4 times the diameter of the observable region; in other words, the 
maximum expected number of copies of the fundamental domain inside the last scatter- 
ing surface in the comoving coordinates is ~8 for compact fiat models without the cosmo- 
logical constant^. Assuming that the primordial power spectrum is scale-invariant(n = 1) 
then the large-angle power is strongly suppressed since fluctuations on scales larger than 
the physical size of the spatial hypersurface are strongly suppressed. 

In contrast, a large amount of large- angle fluctuations can be produced for compact 
low density models owing to the decay of gravitational potential near the present epoch 
which is known as the integrated Sachs- Wolfe effect. If the spatial geometry is flat or 
hyperbolic then the physical distance between two separated points Pi and P2 which sub- 
tends a fixed angle at the observation point O becomes larger as the distance from O 
to Pi or P2 increases. Large-angle fluctuations can be generated at late epoch when the 
fluctuation scale "enters" the topological identification scale L. Recent statistical analyses 
using only the power spectrum have shown that the constraints on the topology are not 
stringent for small CH models including the smallest (Weeks) and the second smallest 
(Thurston) known manifolds and an non-arithmetic orbifold p, ^, 0, 11 1 and also for a 
flat compact toroidal model with the cosmological constant. ||12|| 



These results are clearly at odds with the previous constraints |T^, |TJ] on CH models 
based on pixel-pixel correlation statistics. It has been claimed that the statistical analysis 
using only the power spectrum is not sufficient since it can describe only isotropic (statisti- 
cally spherically symmetric) correlations. This is true inasmuch one considers fluctuations 
observed at a particular point. Because any CH spaces are globally anisotropic, expected 
fluctuations would be statistically globally anisotropic at a particular point. 

However, in order to constrain CH models, it is also necessary to compare the ex- 
pected fluctuation patterns observed at every place of the observer to the data since all 
CH spaces are globally inhomogeneous. |r5[ It should be emphasized that the constraints 



obtained in the previous analyses are only for CH models at a particular observation 



^There is an ambiguity in the definition of L. Here we define L as twice the diameter, i.e. the longest 
geodesic distance between arbitrary two points in the space M . Alternatively, one can define L as twice 
the minimum value of the injectivity radius MiupgM {finjip)}- Injectivity radius rmj at a point p is 
equal to the radius of the largest connected and simply connected ball centered at p which does not cross 
itself. 

^The constraints are for models in which the diameter D is comparable to the injectivity radius rinj ■ 
If D is much longer than rj„j then the constraint on rj„j is less stringent. 
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point P where the injectivity radius is locally maximum for 24 particular orientations. P 
is rather special one in the sense that some of the mode functions often (eigenfunctions 
of the Laplacian) have symmetric structure. [IC] It is often the case that the base point 
P becomes a fixed point of symmetries of the Dirichlet domain or of the space itself. At 
different places the statistically averaged temperature correlations significantly changes 
their anisotropic structure as implied by the random Gaussian behavior of the mode func- 
tions of the Laplace-Beltrami operator. Thus it is of crucial importance to carry out 
Bayesian analyses incorporating the dependence of the statistically averaged correlations 
on the position of the observer. 

On the other hand, recent CMB observations have succeeded in measuring the am- 
plitude and the peak of the angular power spectrum on small angular scales. Assuming 
that initial fluctuations are purely adiabatic seeded by quantum fluctuations then the 
first acoustic peak at I ~ 200 implies nearly flat geometry, [p!^, This seems to deny 
any spherical or hyperbolic models. However, if one considers generalized initial fluctu- 
ations that include isocurvature modes of baryons, cold dark matter and neutrinos, the 
uncertainty in estimating the curvature significantly increases. ||19| In order to decrease the 
uncertainty in the cosmological parameters one needs information of polarization [p!9|, pO 



which will be supplied by the future satellite missions, namely MAP and Planck. Un- 
til then we should interpret the recent CMB observations as which implies "not grossly 
curved spatial geometry" rather than "rigorously flat geometry" . 

Even if the spatial geometry is nearly flat, the effect of the non-trivial topology is still 
prominent provided that the spatial hypersurface is much smaller than the observable 
region at present. However, we cannot expect a very small CH space since there is a 
lower bound Vmin for the volume V of CH spaces. The precise value of Vmin is not known 
but it has been proved that V > 0.1667-R^ where R denotes the curvature radius for CH 
manifolds. The smallest known example has been discovered by Weeks [PT| which has vol- 
ume 0.9427-R^ called the Weeks manifold. If we allow CH orbifold models then the volume 
can be much smaller. For instance, the smallest known orbifold has volume V = 0.0391i?^. 

Suppose that the cosmological constant or the "missing energy" with negative equation- 
of-state dominates the present energy density as the recent observations of SNIa imply 



2^, ^ then the comoving radius of the last scattering surface in unit of curvature radius 
becomes large because of the slow decrease in the cosmic expansion rate in the past. For 
instance, the number N of copies of the fundamental domain inside the last scattering in 
the comoving coordinates is approximately 17.2 for a Weeks model with Q\ = 0.65 and 
f2m= 0.2whereas A^ = 2.5 if r2A = and flm= 0.85. N can be much larger for small orbifold 
models. Thus, even in the case of nearly flat geometry, there are possibilities that we can 
observe the imprint of the non-trivial topology. 

In this paper, we investigate the CMB anisotropy in CH models(manifold and orb- 
ifold) with small volume with or without the cosmological constant. In section 2 we briefly 
describe the necessary ingredients of the geometry and topology of CH spaces. In section 
3 the time evolution of the scalar-type perturbation which can be used for simulating the 
CMB anisotropy is described. In section 4 we estimate the degree of suppression in the 
large-angle power for manifold and orbifold models. In section 5 we carry out Bayesian 
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analyses using the COBE-DMR data. In the last section we summarize our results. 



2 Compact Hyperbolic Spaces 

The discrete subgroup F of the orientation-preserving isometry group of the simply- 
connected hyperbolic 3-space (which is isomorphic to PSL{2, C)) is called the Kleinian 
group. Any orientable CH 3-spaces (either manifold or orbifold) can be described as com- 
pact quotientsQM = H^/T. Let us represent as an upper half space (xi,X2,X3) for 
which the metric is given by 



ds = -2 , (1) 

^3 

where R is the curvature radius. In what follows, R is set to unity without loss of 
generality. If we represent a point p on the upper-half quaternion whose 

fourth component equals zero, then the actions of PSL{2, C) on U C U {oo} take the 
form 

, ap + b 

^■p—,p= ad-bc=l, p = z + X3j, z = xi+X2i, (2) 

cp + d 

where a, b, c and d are complex numbers and 1, i and j are the components of the basis 
of quaternions. The action 7 is explicitly written as 

7 : U C U {00} ^ i7^UCU{oo}, 



f {az + b){cz + d) + acxl X3 \ 

7: ZXi,X2 ,X3 ^ 12 I I 12 2 '1 1 ,12 I I 12 2 h ^ 

\ Icz -I- dp -I- Icl^'xg \cz + d\'^ + \c\^x^ J 

where a bar denotes a complex conjugate. Elements of F for orientable CH manifolds are 
SL{2, C) conjugate to 

which are called loxodromic if 7^ and hyperbolic if = 0. 

There are other classes of isometrics of which have fixed points, namely, parabolic 
and elliptic elements. In matrix representation, they are S'L(2,C) conjugate to 

± f Q 1 ) ' V^'^^^olic (5) 



± ( "^P^^'^/^) _,-^/o^ ) , < < 2vr, elUptzc (6) 



and 

^ ( exp 

exp(i0/2) 

•^The orientation-reversing isometry group of is isomorphic to PSL{2, C)j. If F includes an element 
which is isomorphic to the discrete subgroup of PSL{2, C)j then the compact quotient M is unorientable. 
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, respectively. An orientable CH orbifold can be realized as a compact quotient by F 
which includes an elliptic element. If 7 includes a parabolic element, then the quotient 
space is called a cusped hyperbolic manifold. A cusp point corresponds to a fixed point at 
infinity by the parabolic transformation. In what follows we denote a "cusped manifold" 
as a finite-volume hyperbolic manifold with one cusp point or more. 

Topological construction of CH spaces starts with a cusped manifold Mc- Let us first 
consider the two-dimensional case. If we glue two ideal hyperbolic triangles together 
on each sides, then we have a thrice punctured sphere. Because each side is isometric 
to a line, we can modify the gluing map by an arbitrary translation (in this case it is 
parametrized by M^). However, the obtained surface is not always complete. Let h he a 
horocycleQsegment, that is orthogonal to the two sides of the triangle centered at a vertex 
V. Extend h counterclockwise about v by horocycle segments that meet successive edges 
of ideal triangles. Finally it re-enters the orthogonal triangle as a horocycle concentric 
with h, at a distance d{v). The obtained surface is complete if and only if d{v) = for 
all vertices f l^^. In this case, the condition d{v) = yields a hyperbolic manifold with 
three cusp points where the horocycles correspond to periodic geodesies associated with 
parabolic elements. 

Similarly, one can construct a cusped 3-manifold by gluing ideal polyhedra. Here we 
consider gluing of a set of n ideal tetrahedra Ti,...,T„. Tj is determined by the three 
dihedral angles a, (3 and 7 which satisfy a + /? — 7 = vr. Because the dihedral angles 
of opposite edges are equal, T, is essentially determined by two parameters. In order to 
represent the shape of an ideal tetrahedron Tj , one can use the link L{v) of a vertex v at 
infinity which is an Euclidean triangle realized as the intersection of a horosphere| with Tj 
in the neighborhood of v. L{v) can be described by two parameters and determines Tj up 
to congruence. In order to parameterize L{v), it is convenient to represent the Euclidean 
plane as C.|]2^ To each vertex vi of a triangle (ui, f2, ^3) (the vertices are labeled in a 
clockwise order)in C we associate the ratio 

f2 - Vi 

of the sides adjacent to vi. The ratios(with any starting point) should satisfy 

I - zivi) + z{vi)z{v2) = 0. (8) 

Thus any z{vi) determines the other two z{vj)s. An ideal tetrahedron can be specified by 
any z{vi) having two parameters up to orientation-preserving similarity (figure |I|). The 
gluing conditions of ideal tetrahedra Ti, ...,T„ are algebraically given by 

z{ei) ■ z{e2) z{em) = 1, 

argz(ei) H h arg2;(em,) = 27r, < arg < vr (9) 



^Horocycles centered at a vertex v at infinity are the curves orthogonal to all lines through v. 
^Horospheres centered at a vertex v at infinity are the surfaces orthogonal to all lines through v. 
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Figure 1: Parameterization of an ideal tetrahedron. Each edge e is labelled with a com- 
plex number z{e) and opposite edges have the same label. The four links (=Euclidean 
triangles realized as the intersections with horospheres in the neighborhood of vertices) 
are congruent each other. An ideal tetrahedron up to orientation-preserving similarity 
can be specified by one link which is parametrized by one of the z{e). 



where e^s are the edges of ideal tetrahedra. Note that even if the above conditions are 
satisfied the obtained manifold is sometimes incomplete where the developing image of 
L{v) is C — after an appropriate translation. If L{v) properly tessellates C then the 
obtained manifold M is complete and the boundary of the neighborhood of a vertex at 
infinity is either a torus or a Klein bottle. M = is topologically equivalent to the 
complement of a knot K or link L (which consists of knots) in a 3-sphere or some other 
closed 3-space. A surgery in which one removes the tubular neighborhood N of K whose 
boundary is homeomorphic to a torus, and replace by a solid torus so that a meridian^ 
in the solid torus goes to {p, q) curve]] on A^ is called {p, q) Dehn surgery. Except for a 
finite number of cases, Dehn surgeries on K where p and q are co-prime always yield CH 
3-manifolds which implies that most compact 3- manifolds are hyperbolic [p5|] since every 



closed 3-manifolds can be obtained by performing such Dehn surgeries. CH 3-orbifolds 
whose singular sets consist of lines can be obtained by {p, q) Dehn surgeries where p and 
q are not co-prime. 



^Given a set of generators a and b for the fundamental group of a torus, a closed curve C which 
connects a point x in the torus with ax is called a meridian and another curve which connects a point x 
with bx is called a longitude. 

^If a closed curve C connects a point x with another point {pa + qh)x where p and q are integer, C is 
called a (p, q) curve. 
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Figure 2: The Dirichlet fundamental domains of the Weeks manifold (volume= 
0.9427) (left) and the smallest orbifold that is covered by the Weeks manifold (volume= 
0.0786) (right) viewed from opposite directions in the Klein(projective) coordinates. Col- 
ors on the faces represent the identification maps. 



The above procedure for obtaining CH 3-spaces is manifestly defined but somewhat 
complicated for computation by hand. A computer program "SnapPea" developed by 
Weeks [ p^] can numerically handle the procedure, namely computation of hyperbolic struc- 
ture of cusped manifolds, performing various Dehn surgeries and so on. SnapPea can also 
compute the volume, fundamental group, homology, symmetry group, Chern-Simon in- 
variant, length spectrum and other quantities. 

In contrast to compact fiat spaces with no scale, there is a lower bound for the vol- 
ume of CH 3-spaces. It has been proved that the volumes of CH 3-manifolds cannot be 
smaller than 0.16668. . .|2^ although no concrete example with such volume is known. 
The smallest known CH 3-manifold is called the Weeks manifold with volume 0.9427 
which can be obtained by a (3,-1) Dehn surgery on a cusped manifold m003 with volume 
2.0299(one of the smallest known cusped manifold). Performing a (-2,3) Dehn surgery 
on m003 yields the second smallest known manifold called the Thurston manifold with 
volume 0.9814. Except for a finite number of cases, Dehn surgeries on m003 yield CH 
manifolds with volume less than that of m003.[|25[ As \p\ and |g| becomes large, the vol- 
ume of the space converges to that of the original cusped manifold Mc- Therefore, for 
the purpose of classification of CH spaces it is better to exclude "almost non-compact" 
spaces that are very similar to Mc.| The volume of CH 3-orbifolds can be much smaller 
since every CH 3-orbifold is covered by a CH 3-manifold (there are an infinite number 
of coverings). For instance, the smallest volume of a CH 3-orbifold whose finite-sheeted 
cover is the Weeks manifold is V = 0.9427/12 = 0.0786 (see figure 0)since the isome- 
try group of the Weeks manifold is D6 (dihedral group), which does not have any fixed 



*The Hodgson- Weeks census consists of the data of 11,031 CH manifolds with length of the shortest 
periodic geodesic I > 0.3 and volume V < 6.4535. 
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planes. The smallest known orientable 3-orbifold has volume V = 0.0391 which also does 
not have any fixed planes. ||2^ A known lower bound for the volume of orientable orbifolds 
isV > 0.0000013. 



3 CMB Anisotropy 

In what follows we assume that the matter consists of two components:a relativistic and 
an non-relativistic one {Qm = + ^n)- From the Friedmann equation, the integral 
representation for the time evolution of the scale factor a (normalized to 1 at present 
time) in terms of the conformal time r] is given by 

r]{a) = R I K^O, (10) 

Jo VsignK((fi„a + + nAa'^)/{^m + ^^a - 1) - a^) 

where R = \K\~^/'^ is the curvature radius at present. The above presentations can also 
be written in terms of some special functions but they are too complex to describe here. 

The scalar perturbation equations are exactly same as those without the cosmological 
constant. Let us consider the adiabatic case and assume that the anisotropic pressure 
is negligible. Then from the linear perturbation equation, the time evolution of the 
Newtonian curvature perturbation $ of mode is given by, |^ 

+ 3n{i + c2)$' + c2fc2$ + {2n' + (1 + 3cl){n^ - k))^ = o, (ii) 

where ' denotes the conformal time derivative d/drj, Ti = a'/a, G is the Newton's constant, 
Cs is the sound speed of the fiuid which is explicitly written in terms of the density of 
relativistic matter and non-relativistic matter 

cl = (12) 

O Pr -r ^pm 

Let us consider the early universe when the radiation is dominant and the effect of the 
curvature and the cosmological constant is negligible. Then ( p!OD yield the evolution of 



the scale factor as a oc and (|TT]) is reduced to 



4 fc2 

$" + + — $ = 0, (13) 
•q 3 

which has an analytic solution 

$(//) = Tj^^ {{ijjTj cos{ujvi) — sm{ijjrii))Ci + {ujrj sm{ijjvi) — cos{ujvi))C2}, (14) 
where uj = k/ \/3 and Ci and C2 are constants which depend on k.'^^ Around r/ = 0, (|1^) 



can be expanded as 

$(r^) = r]-\C2 + - ^t^V - + 0{{uoT]f)}. (15) 



^2 22 ^1 , ,3^3 ^2 4 4 
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Thus for the long- wavelength modes(a;?7 << 1), the non-decaying mode (C2 = 0) has 
constant amplitude and the time derivative at the initial time vanishes, 

$ ~ const, $'(0) = 0. (16) 

During the radiation- matter equality time the equation of state change from w = 1/3 to 
w = 0. Then the amplitude of long-wavelength modes changes by a factor of 9/10. In the 
curvature or A dominant epoch the amplitude decays as 1/a. 

Assuming that the matter is dominant at the last scattering ^ and the anisotropic 
pressure is negligible, the temperature fluctuation on large angular scales can be written 
in terms of the Newtonian curvature perturbation $ as 

^(„, . (,„ _ _ 2 £" ^^'■'•'g;-')") ... (17) 

where n denotes the unit vector which points the sky and rj^ and tjq correspond to the 
last scattering conformal time and the present conformal time, respectively. [^, |3^ The 
first term in the right-hand side in (0) describes the ordinary Sachs- Wolfe (OSW) effect 
while the second term describes the integrated Sachs- Wolfe (ISW) effect. 



4 Angular power spectra 

In order to analyze the temperature anisotropy in the sky it is convenient to expand it in 
terms of spherical harmonics Yim as 

AT 

^(n) = ^aimYimin), (18) 

Im 

where n is the unit vector along the line of sight. Because each mode of scalar perturbation 
evolves independently in the locally isotropic and homogeneous background space, we 
need only the information of eigenmodes Ui, of the Laplace-Beltrami operator A in the 
CH space (i.e. elements of an orthonormal basis of LF'{H^ /V)) which can be expanded 
in terms of eigenmodes in the simply connected infinite hyperbolic space in the spherical 
coordinates (x, 0, (p) as 



Im 



r(/ + 1 + m) / I ^p-/-i/2, 
[ui) Y sinhx 



= ^-F7-^A/-rT-^.ri/2 (coshx), (19) 
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For low density models, the ordinary Sachs- Wolfe term in (h7h may not be valid since the last 



scattering may occur before full-matter domination. However, on large angular scales (17) still gives a 
good approximation since the integrated Sachs- Wolfe term dominates over the ordinary Sachs- Wolfe term 
in low density models. 
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where = fc^ — 1, Yim denotes (complex) spherical harmonics, P is the associated 
Legendre function and ^i,im represents the expansion coefficients. 

Assuming that the initial fluctuations obey the Gaussian statistic, and neglecting the 
tensor-type perturbations and anisotropic pressure and matter is dominant at the last 
scattering epoch, the angular power spectrum Ci =< la^mP > (<> denotes an ensemble 
average over the initial perturbation) for CH models can be written as 



1 r° d^t 

FuiiVo) = -Mv*)Xui{Vo-V*)+yd7]^X,i{r]o-r]), (20) 



where u = a/A?— T, ^$(1^) is the initial power spectrum, and rj^, and rjo are the last 
scattering and the present conformal time, respectively. $j = $/$(0) and V{M) is the 
comoving volume of the space. The (extended) Harrison-Zel'dovich spectrum corresponds 
to V.^iy) = Const, which we shall use as the initial condition. It should be noted that 
C;'s depend on the position and orientation of the observer since CH spaces are globally 
inhomogeneous. Therefore, it is better to consider the ensemble average Ci =< Ci > 
taken over the position of the observer. Using Weyl's asymptotic formula one can easily 
show that C;'s coincide with those for the infinite counterpart in the limit u ^ 00 provided 
that < l^j/ZmP >oc Q Thus all we have to do is to numerically compute the expansion 
coefficients ^^im which contain the information of global topology and geometry.|^ 

In order to estimate the suppression in the angular power Ci it is convenient to write 



Ci in terms of the transfer function Ti{k) which satisfies 



^Ic, = / r,=(*)Pwf , (21) 

where Ci denotes the angular power for the infinite counterpart. For CH models, Ci is 
written as a sum of Tf{k) for discrete wave numbers fcj. 

On large angular scales, the behavior of the angular power Ci is determined by low- 
lying modes which are susceptible to the global topology of the background geometry 
since the amplitude of each mode is proportional to ~ z/^^l-Fj^/p. As shown in figure ^ 
one can see that the transfer function corresponding to the first eigenmode T^{ki) mimics 
the angular dependence of Ci. 

The large-angle power owing to the OSW effect suffers a significant suppression since 
fiuctuations on scales beyond the actual size of the space are not allowed. The angular 
cutoff scale Icut above which the OSW contribution suffers a suppression corresponds to 

Although we have assumed that the normalization factor in < |<^iy;mP > does not depend on the 
volume of the space, the volume factor in (20) might not be necessary if < |^iy;mP > is proportional to 



l/V{M). 

^^The low- lying eigenmodes of the two smallest manifolds (Weeks and Thurston) have been computed 
using the direct boundary element method. ||3^ For other small CH manifolds the eigenvalues have been 
computed using the periodic orbit sum method. [p4[ 
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the angular scale of the first eigenmode on the last scattering surface, which can be written 
in terms of the comoving volume V{M) and the smallest non-zero wavenumber ki as 

lent = J Vr-i((sinh(2(i?i55 + ra.e)) ' smh{2{RLss - rave)) - 4r,,e)), (22) 

where V = 7r(sinh(2rat,e) — 2rave) and Rlss denotes the comoving radius of the last 
scattering surface. For the Weeks models {V = 0.94 and ki = 5.27), Icut = 26 for 
Qq = 0.2 and Icut = 7 for fio = 0.6. From figure ^ one can see that the angular scales Icut 
corresponding to the intersection of k = ki and the OSW ridge of the transfer function 
well agree with these analytic estimates. On smaller angular scales / >> Icut, the powers 
asymptotically converge to those for the infinite counterpart. 

Because the large-angle power in the COBE-DMR data is nearly flat, it seems that the 
suppression leads to a bad fit to the data on large angular scales. However, for low density 
models we should consider the ISW effect owing to the gravitational potential decay at the 
curvature dominant epoch 1 + 2; ~ (1 — no)/f2o or A dominant epoch 1 + 2; ~ (f^A/f^o)^^^ 
well after the last scattering time. For a given angular scale / the comoving scale k~^ of 
a fluctuation that is produced at late time is decreased. On the other hand, fluctuations 
do not suffer suppression if the comoving scales of fluctuations are sufficiently smaller 
than the actual size of the space k^^ « Tave- Therefore, the suppression on the ISW 
contribution is less stringent compared with the OSW contribution. Interestingly, the 
suppression on the angular power owing to the mode-cutoff reduces the excess power 
owing to the ISW effect, resulting in a nearly flat power with a slight suppression on large 
angular scales 2 < / < 10. 

The suppression of the power is crudely estimated by the number of the copies 
of the fundamental domain inside the last scattering surface in comoving coordinates. 
Suppose that f2A = 0. Then the comoving radius of the last scattering surface in terms of 
the present curvature radius R 

i?L55~ i?cosh-^(2M„-l), (23) 

gives the comoving volume vlss of the ball inside the last scattering surface 

VLSS ~ TiR\sm\i{2RLss/R) - 2Rlss/R). (24) 

For example, v^ss ~ 490i?'^ for {Qm,^A) = (0.2,0) whereas vlss ~ for (^^^^a) = 
(0.9,0). Thus in nearly flat models, the imprint of the non-trivial topology is prominent 
only for the case where the volume is smaller than R^. However, if one includes the 
cosmological constant then Rlss (in unit of R) becomes large because of a slow increase 
in the cosmic expansion rate in the past. For instance, = 8.7 for a Weeks model (the 
smallest known manifold) with Q\ = 0.7 and Qm =0.2 whereas A^ = 1.2 for Q\ = and 
Qm = 0.9. Furthermore, if one allows orbifold models, A^ can be much larger than these 
values. For instance, A^ = 51.8 for the smallest orbifold (volume= 0.0391i?^) with Q\ = 0.7 
and flm = 0.25 

Now let us consider the angular power spectrum for orbifold models. Assuming no 
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Figure 3: Suppression in the large-angle power AT^ = l{l + l)C'//(27r) for 5 CH mod- 
els (name, volume) = A:(m003(3,-l)(Weeks),0.94), B:(m010(-l,3),1.9), C:(m082(-2,3),2.9), 
D:(m288(-5,1),3.9) E:(s873(-4,1), 4.9) in comparison with the one for the infinite hyper- 
bolic model(denoted as O). All the plotted values are normalized by AT20 for the infinite 
hyperbolic model. The I dependence of AT; can be approximately given by the transfer 
function lT^{ki). First eigenvalues ki for the five models are represented as horizontal 
lines in lower figures. The unit of k is equal to the inverse of the present curvature radius. 
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Figure 4: Suppression in large-angle power AT; = + l)C*;/27r for small (hypothetical) 
orbifold models with = 0.2 and fl\ = 0.7 in comparison with the infinite hyperbolic 
model with the same density parameters. The amplitude is normalized by AT20 for the 
infinite hyperbolic model. 



very short periodic geodesies and no supercurvature modes, one can crudely estimate 



the low-lying eigenvalues |^ from the asymptotic form of the number function(spectral 
staircase) 

n{E) = ^^u^ + C2U^ + cii/ + Co + 0(6-"^^/^), u^ = E-l, (25) 

where n{E) denotes the number of eigenvalues of the Laplace-Beltrami operator equal to 
or less than E and q's are the constants determined from the plane refiection, elliptic and 
inverse elliptic elements of the discrete isometry group and the area of the fixed planes. 
For orientable orbifolds with only fixed lines, the constants vanish except for ci which is 
written in terms of elliptic elements. If we consider only such orbifolds with small volume 
(large Vi) then the dominant contribution comes from the first term in the right hand 
side in (PB]). If we further assume that the eigenmodes have the same pseudo-random 
property as those of CH manifolds then we can readily estimate the angular powers. 



As shown in figure ^ the suppression of the large-angle power for small orbifold models is 
still prominent in the case of nearly fiat geometry. 

Another feature of the non-trivial topology is an increase in the cosmic variance which 
is attributed to the global inhomogeneity and anisotropy of the background geometry. If 
we assume the pseudo-random Gaussianity of the expansion coefficients ^yim as observed 
in the smallest CH manifolds p3|, |T6| then uncertainty in the power (which may be called 
the geometric variance) can be easily estimated. For CH models, aim is given by a sum 



13 



5 10 15 20 

I 



Figure 5: Fractional uncertainty in the angular power ACi = (Var(C*/))^/^/ < Ci > 
for the Weeks models relative to the values for the Gaussian models AC; (Gauss) = 
(Var{Ci)y/y < Ci >= (2/2/ + 1)1/2 

of products of the initial perturbation $^(0) times the expansion coefficients 

aim = 5^$.(0)e.zmi^.i. (26) 

If we fix the values of the primordial fluctuations then aim^s behave as if they are random 
Gaussian numbers and the angular power obeys the distribution. Because the degree 
of freedom {k, m) in the sum is always larger than 2/ + 1, the geometric variance should be 
always smaller than the "initial variance" owing to the uncertainty in the initial conditions. 
Therefore, we expect that the net cosmic variance (aiinitial variance + geometric variance) 
is not significantly greater than the values for the infinite counterpart. 

In order to confirm the validity of the Gaussian assumption for computing the cosmic 
variance, firstly, we have compared the fractional uncertainty in the angular power (2 < 
/ < 20) of the Weeks models using only the lowest 33 numerically computed eigenmodes 
{k < 13) to those using the Gaussian random approximation for the expansion coefficients 
(corresponding to the same eigenvalues). It has turned out that the errors lie within several 
per cent for models with Vim = 0.2,0.4 and 0.6(the relative errors are 0.08 for I = 2 and 
< 0.04 for / > 2). 

Next, taking the contribution of higher modes into account, we have computed the 
fractional uncertainty in the angular power, which can be estimated by using the Gaussian 
approximation for the expansion coefficients (corresponding to the computed eigenvalues 
for A; < 13 and the approximated eigenvalues obtained from Weyl's asymptotic formula 
for k < 60) for the same parameters and we have compared the values with those using 
only the lowest 33 eigenmodes {k < 13) and it has found that the errors lie within several 
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per cent (however, for high matter density models the systematic increase in the variance 
owing to the artificial mode cut-off is much significant as the number of modes that 
contribute to the sum grows). As shown in figure ^, on large angular scales the fractional 
uncertainty in the angular power (for the Weeks models) with low matter density increases 
just 5 to 20 percent relative to the one for the Gaussian models. Thus the previous claim 
that the C^'s have large cosmic variances ||14|| is not correct. On small angular scales, an 
increase in the number of eigenmodes that contributes to the power leads to a decrease 
in the geometric variance. Consequently, the cosmic variance converges to the one for the 
Gaussian model as implied by the central limit theorem. 



5 Bayesian analysis 

In this section, we study the likelihoods of CH models using the COBE data. Although 
here we only study manifold models, we expect the result for orbifold models (with the 
same volume) will not grossly change since the statistical property of eigenmodes are 
expected to be similar with those of manifolds |10|. 

The covariance in the temperature at pixel i and pixel j in the sky map is given by 

Mij =< T.{Tj >= 5^ < aimavm' > WiW;Yim{ni)Yi,^'{nj)+ < N,Nj >, (27) 
I 

where <> denotes an ensemble average taken over all initial conditions, positions and 
orientations of the observer]^ Tj represents the temperature in pixel i, is the experi- 
mental window function that includes effects of beam-smoothing and finite pixel size, Ui 
is the unit vector towards the center of pixel i and < NiNj > is the noise covariance 
between pixel i and pixel j. If the fiuctuations in the sky form an isotropic Gaussian field 
then the covariance is written as 

Mij = i- ^(2/ + l)W^CiPi{ni ■ nj)+ < N,Nj >, (28) 
I 

where Pi is the Legendre function. Assuming a uniform prior distribution for a set of 
cosmological parameters, the probability distribution function of a power spectrum Ci is 
given by 

A(Ci\f) oc ^ exp(-f^-M'\Ci)-f\ (29) 

^ ' ' det^/'M(Q) "^{2 y ^) )^ y > 

where T denotes an array of the data of the temperature at pixels. 

In the following analysis, we use the inverse-noise-variance-weighted average map of 
the 53A,53B,90A and 90B COBE-DMR channels. To remove the emission from the galac- 
tic plane, we use the extended galactic cut (in galactic coordinates).]^ After the galactic 
cut, best-fit monopole and dipole are removed using the least-square method. To achieve 

^^Here we assume that we do not know anything about the position and orientation of the observer. 
The covariance is defined for an isotropic and homogeneous ensemble of observers. 
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Figure 6: Plots of the ratio of likelihoods (assuming that the initial condition is fixed) 
for two smallest CH models (Weeks and Thurston) to a likelihood for the Einstein-de 
Sitter model {Qm = 1-0) with scale invariant initial spectrum {n = l). All the likelihoods 
are marginalized over the normalization of the power. Here we have assumed $^(0) oc 
l/(z/(z/2 + 1)) for CH models. 



efficient analysis in computation, we further compress the data at "resolution 6" (2.6°)^ 
pixels into one at "resolution 5" (5.2°)^ pixels for which there are 1536 pixels in the ce- 
lestial sphere and 924 pixels surviving the extended galactic cut. The window function is 
given by Wi = GiFi where Fi are the Legendre coefficients for the DMR beam pattern p9| 
and Gi are the Legendre coefficients for a circular top-hat function with area equal to 
the pixel area which account for the pixel smoothing effect (the effect of the finite pixel 
size is non-negligible since the COBE-DMR beam FWHM is comparable to the size of 
"resolution 5" pixels). To account for the fact that we do not have useful information 
about the monopole and dipole anisotropy, we set Cq = Ci = lOOmK"^ which renders the 
likelihood insensitive to monopole and dipole moments of several mK. We also assume 



that the noise in the pixels is uncorrelated which is found to be a good approximation. pT 



First of all, we set the initial condition as $^(0) oc l/(z/(z/^ -|- 1)) in order to ap- 
proximately estimate the effect of the suppression in the large-angle power owing to the 
non-trivial topology (here we do not consider $i/(0) as random numbers). Then the 
fluctuations form a pseudo-Gaussian random fleld assuming pseudo-Gaussianity of the 
eigenmodes. As shown in flgure ^ for a wide range of parameters {flm + > 0.1) 
the likelihoods for the smallest CH manifold models (Weeks and Thurston) are better 
than one for the Einstein-de Sitter model with scale-invariant spectrum [n = 1) where 
AT; = (/(/ + l)Q/27r)V2 is almost constant in /. One can see the better flts to the COBE 
data for three parameter regions:l. Qm = 0.5 ~ 0.7 with small for which the angular 
power is peaked at / ~ 4 which corresponds to the flrst ISW ridge of the transfer function; 
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Figure 7: Plots of angular power spectrum AT; = (/(/ + l)Ci/2nY/'^ of the COBE-DMR 



data (box) |^ in comparison with those for the Thurston models with ilm = 0.2, 0.6 and 
a flat A model with Qm = 0.3 and Q\ = 0.7. 



2. ^Irn = 0.85 ~ 0.9 where the angular scale I ~ 4 which corresponds to the SW ridge;3. 
Qrn ~ 0.2 where the slope of the power on large angular scales / < 10 fits well with the 
data. 

For low matter density models with infinite volume (that is simply connected) the ISW 
effect leads to an excess power on large angular scales. Therefore the fit to the COBE 
data is not good because of the low quadrapole moment in the data(figure |^. In contrast, 
for small CH models, as we have seen, the excess power owing to the ISW effect is miti- 
gated by a suppression owing to the mode-cutoff of the eigenmodes. Therefore, hkelihoods 
for small CH models with low matter density are significantly improved compared with 
the infinite counterparts (figure |^). As the volume increases the likelihoods converge to 
those of the infinite counterparts although the convergence rate depends on cosmological 
parameters. One can see in figure ^ that the conspicuous difference for Q^n = 0.1 — 0.3 
still persists for volume~ 6 whereas such difference is not observed for nearly fiat cases 
{Qtot = + = 0.9). Roughly speaking, the difference in the power depends on the 
number of the copies of the fundamental domain inside the observable region at present. 

Next, we consider the effect of the non-diagonal elements which we have neglected so 
far. The likelihood for a homogeneous and isotropic ensemble is obtained by marginalizing 
the likelihoods all over the positions Xobs and the orientations a of the observer, 

^= y/g dyiobsdaX{y.obs, a), (30) 
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Figure 8: Plots of the ratio of likelihoods for the infinite hyperbolic "open" models {n = 1) 
to one for the Einstein-de Sitter model (n = 1) (left) and likelihoods for the Thurston mod- 
els relative to one for infinite hyperbolic models (right). All the likelihoods are marginal- 
ized over the normalization of the power. Here we have assumed $^(0) oc l/(z/(z/^ + 1)) 
for the Thurston models. The slight improvement in the likelihood of infinite hyperbolic 
models with Qm < 0.1 is caused by the absence of the supercurvature modes. For the 
Thurston models with Qm = 0.1 ~ 0.3 the likelihoods are significantly improved. 



where we assume a constant distribution for the volume elements ^(hi.ohs and da of a Lie 
group 50(3) with a Haar measure. Assuming that the initial fluctuations are Gaussian 
with Harrison-Zel'dovich spectrum {n = 1), the likelihood A(xo;,<., a)) is given by (123) and 



( p9D where < aimO-i'm' > is written as 



m' FylPvV 1 ('^-^) 

V^W^ + 1) 

where <> denotes an ensemble average taken over the initial condition and Fyi describes 
contribution from the OSW effect and the ISW effect, respectively. Note that ^^im's are 
functions of :s.obs and a. In order to compute likelihoods, we use a compressed data at 
"resolution 3" (20.4°)^ pixels in galactic coordinates for which there are 60 pixels surviving 
the extended galactic cut for efficient analysis in computation. Although the information 
of fluctuations on small angular scales / > 10 is lost, we expect that they still provide 
us sufficient information for discriminating the effect of the non-trivial topology which is 
manifest on large-angular scales. 

For the Thurston model, it has turned out that only 0.09 percent of the total of 
1500 positions with 10 orientations are larger than the mean value. We have also com- 
puted likelihoods for 5100 positions with 40 orientations. Then the percentage has re- 
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Figure 9: Plots of the ratio of approximated likelihoods for 27 CH models with volume 
(0.94 — 6.15) to a likelihood for the infinite hyperbolic counterparts. All the likelihoods 
are marginalized over the normalization of the power. Here we have assumed $^(0) oc 
l/(i/(z/2 + 1)) for CH models. 

duced to 0.02 from 0.09. The ratio of the likelihood marginalized over 5100 positions 
and 40 orientations to the likelihood of the infinite hyperbolic model with Qm = 0.4 is 
Logj^Q (A/A (open)) = —1.4 and the maximum value is Log]^Q(A(max)/A(open)) = 3.6. For 
17 cases out of 204000 realizations, the likelihoods are much better than one for the infinite 
counterpart. Although 99.98 percent choices of position and orientation are ruled out, the 
likelihood for the remaining choices is approximately 4000 times larger than that for the 
infinite counterpart which boosts the probability of having skymaps consistent with the 
data. It has turned out that the positions that give a better fit to the data are scattered 
in the manifold and do not coincide with the point where the injectivity radius is locally 
maximal(=the center of the Dirichlet domain). This suggests that we are accidentally 
put at a certain place with a certain orientation. In other words, the observed sky gives 
us partial information about the position and orientation of the observer in the mani- 
fold although it is not enough to determine the values uniquely. The best-fit quadrupole 
normalization is Q = (5C2/47r)^/^ = 0.022 ± 0.0024:TnK for 50 choices of position and 
orientation of the observer that satisfy A/A(open) > 0.03 while Q = 0.027 ± 0.0052mK 
for a total of 15000 realizations. For "bad" choices the best-fit normalization is somewhat 
high since it gives a large cosmic variance whereas the normalization is much lowered for 
"good" choices. One can see in figure ^ that a random realization of the skymap for 
the Thurston models does not appear grossly inconsistent with the COBE-DMR data. It 
turns out that the statistically averaged anisotropic correlation pattern depends sensitively 
on the position of the observer. 
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Figure 10: Plots of the ratio of "rigorous" likelihoods incorporating the effect of inhomo- 
geneity and anisotropy of the background geometry for two compact multiply connected 
models relative to one for the infinite counterparts in ascending order(full curve) for 2000 
random realizations of orientation for a closed fiat toroidal (where the Dirichlet fundamen- 
tal domain is a cube with sides L = 2.0Hq^) model with Q^n = 0.1, f^A = 0.9 (top) and 
for 1500 random realizations of position and 10 realizations of orientation for a Thurston 
model with flm = 0.4 (bottom). The dashed lines denote the ensemble averaged values. 
The number of the copies of fundamental domain inside the observable region is 64.7 
for the compact flat toroidal model and 72.3 for the Thurston model. 
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Figure 11: The top map shows the inverse-noise- variance-weighted averaged map of the 
COBE-DMR data at 53GHz and 90GHz after Wiener filtering assuming an angular power 
spectrum (averaged over the position of the observer) for a Thurston model with = 0.4. 
Each of the mid and the bottom map shows a random realization of the CMB anisotropy 
(convolved with the COBE-DMR beam) at a position that gives a best-fit to the COBE- 
DMR data. For all the maps the "extended galactic cut" proposed by the COBEl-DMR 
team was used and both the monopole and dipole were removed. 



21 



The likelihood analyses in ||13| , are based on correlations for a particular choice of 
position P where the injectivity radius is locally maximal with 24 orientations. It should 
be emphasized that the choice of P as an observing point is very special one. For instance 
it is the center of tetrahedral and Z2 symmetry of the Dirichlet domain of the Thurston 
manifold although they are not belonging to the symmetry of the manifold. [0 In mathe- 
matical literature it is standard to choose P as the base point which belongs to a "thick" 
part of the manifold since one can expect many symmetries. However, considering such a 
special point as the place of the observer cannot be verified since it is inconsistent with the 
Copernican principle. Because any CH models are globally inhomogeneous, one should 
compare fluctuation patterns expected at every place with every orientation. 

To see the influence of global inhomogeneity, it is illustrative to compare the relative 
likelihoods(to the infinite counterpart) of the Thurston model with those of compact fiat 
toroidal model (obtained by gluing opposite faces of a cube by three translations) which 
is globally homogeneous (i.e. every action of the discrete isometry group is a Clifford 
transformation). As shown in figure 10, for the toroidal model(that has approximately 
the same proportion to the currently observable region in size in comparison with the 
Thurston model), dependency on choices of position of the observer is less significant. 
Therefore one does not need a number of realizations for estimating the likelihood. In 
contrast, for CH models, one needs a sufficient number of realizations since the proportion 
of choices of position and orientation of the observer that give a better fit to the data is 
considerably small. Thus if one treats CH models like the toroidal model [l^, one 
gets misleadingly small values of the likelihoods. 

Although our result is based on the numerical computation, it is a natural result if one 
knows the pseudo-random behavior of eigenmodes of CH spaces. For each choice of posi- 
tion and orientation of the observer, a set of expansion coefficients ^yim of eigenmodes is 
uniquely determined (except for the phase factor), which corresponds to a "realization" of 
independent random Gaussian numbers. By taking an average over the position and the 
orientation, the non-diagonal terms proportional to < iuimiuVm' >,l ^ V ,m ^ m' vanish. 
In other words, a set of anisotropic patterns all over the place in a CH space comprises 
an almost isotropic random field. Consider two realizations A and B of such an isotropic 
random field. The chance you would get an almost similar fiuctuation pattern for A and 
B would be very low but we do have such an occasion. Similarly, the likelihood at a 
particular position with a certain orientation is usually very low but there are cases for 
which the likelihoods are considerably high. Thus we conclude that the COBE constraints 
on small CH models are less stringent as long as the Gaussian pseudo-randomness of the 
eigenmodes holds. 



6 Summary 

In this paper, we have explored the CMB anisotropy in small CH models with or without 
the cosmological constant. Assuming adiabatic initial perturbation with scale- invariant 
spectrum (ra = 1), the angular power owing to the OSW effect suffers a prominent suppres- 
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sion since fluctuations beyond the size of tlie fundamental domain at the last scattering 
are strongly suppressed. However, for low matter density models, the suppression in the 
large-angle power is less stringent because of the significant contribution from the ISW 
effect caused by the decay of the gravitational potential at the A or curvature dominant 
epoch. A slight suppression in the large-angle power in such models explains rather natu- 
rally the observed anomalously low quadrupole which is incompatible with the prediction 
of the standard Friedmann-Robertson- Walker models. 

As we have seen, the likelihood of CH models (assuming pseudo-Gaussianity of eigen- 
modes) depends sensitively on the choice of orientation and position of the observer. 
Because the likelihood marginalized over the orientation and position (assuming equal 
probability for each choice) is comparable to the value of the infinite counterpart, we 
conclude that constraints on CH models are less stringent. It should be emphasized that 
the dependence of the likelihood on the position of the observer is of crucial importance 
which has been ignored by previous literature. Closed multiply connected constantly 
curved 3-spaces that are globally homogeneous are limited to some spherical spaces and 
flat 3-tori.|T^ For "bad" choices of the position and orientation, the best-fit amplitude 
tends to have a large value (allowing a large variance), leading to misleadingly strin- 
gent constraints on the models. Surprisingly, the statistically averaged anisotropy of the 
correlation seems to disappear if marginalized all over the place which is related to the 
pseudo-random property of eigenmodes. 

Even in the case of nearly flat geometry, the signature of the non-trivial topology is 
still prominent if the space is sufficiently small compared with the observable region at 
present. If we allow orbifold models then the volume can be much smaller than manifold 
models. We have seen that a slight suppression in the large-angle power is still prominent 
for orbifold models with volume 0.01 — O.IR^ for Qm = 0.2 and f^A = 0.7 which is con- 
sistent with the result in. [|rT|] However, the existence of singularities (where the curvature 
diverges) may cause some problems for any orbifold models. If plane-like singularities 
were present (e.g. tetrahedral orbifolds) astronomical objects with peculiar velocity would 
easily collide with the plane (we may call such a model as a "billiard universe"). On the 
other hand, the observational effects caused by the presence of fixed lines or "strings' 
might be less prominent and much safer. 
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